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ABSTRACT 


This report presents an approach for improving the stability 
of load-frequency for single area and multi-area interconnected 
power systems using P.I.D. controller i.e. proportional, integral 
and derivative controls. These three components of the P.I.D. 
controller have been used in various combination in forward and 
feedback paths. Only some of the combinations were found to give 
improved transient and steady state responses. The proposed 
technique is applied to single area power system as well as 
two-area interconnected power systems like thermal-thermial and 
hydro-thermal systems. The P.I.D. controller can be applied to 
multi-area-interconnected power systems. 

A comparison of the Load Frequency Control (LFC) response 
with Conventional and Variable-Structure System (VSS) control and 
proposed P.I.D. control strategies shows that, with the 
application of the proposed algorithm, the system performance is 
improved significantly. 
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CHAPTER - 1 


IHTRODUCTI<»l 


1.1 CENERAL 

It is important to maintain the frequency constant in the 
power systems. Variation of the frequency in the power 
system may cause problems. Some of these problems are listed 
be lows 

(i) Most AC motors run at speeds that are directly 
related to the frequency. The speed and induced electromotive 
force (EMF) may vary because of the change of frequency of the 
power circuit. 

(ii) When operating at frequencies below 49.5Hz, some types 
of steam turbine, certain rotor state undergo excessive 
vibrations which results in metal fatigue and blade failures. 

(iii) The turbine regulating devices fully open when the 
frequency falls below 49H2 . This situation causes extra loading 
on the generators. A decrease in frequency may cause reduction 
in the equipment efficiency, e.g., feed pumps. 

(iv> A large number of electrically operated clocks are 
used fof certain purposes. If frequency changes, it can 
influence the operation of these clocks. 
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(v) The change in frequency can cause fnaloperalion of 
power convertors by producing harmonics. 

The load-f requecy control is installed for better exchange of 
Mega-Watt Hour <M. W.H) , stabi 1 i ty and fTiinimize above mentioned 
problems. For reasons of economy and system reliability, 
neighbouring power systems are interconnected, forming an 
augmented system refered to as ’power pool’. The net power flow 
on the tie-lines connecting a system to the external system is 
frequently scheduled by apriori contact basis. System 
disturbances caused by load fluctuations result in changes in 
tie-line power and system frequency which gives rise to a "load 
frequency control" (LFC) problem- 

To study the load frequency control, the excitation control 
influence on the system performance can be ignored because of the 
following reasonss 

(i) For small changes, it can safely be assumed that 
change in real power only causes change in load angle- Hence, 
the variations in load angle causes momentry change in generation 
speed. This means generator has been supplied with sufficient 
reactive power to maintain the voltage constant at the terminals- 

(ii) The major time constant involved in the excitation 
system is generator’s field time constant. The time constants 
involved in load frequency control are those of the turbine and 


inertia of generator. The Latter time constants are much longer 

compared to the excitation system time constant. Hence it is 

fair to assume that the transients of the excitation system will 
vanish much faster then the transients of the load frequency 
control system and does not affect the response of the load 
frequency control system. 

The load frequency control is based on an error signal 
called Area Control Error (ACE) which is a linear comnbination of 
net interchange and frequency errors. The conventional control 
strategy used in industries is to take the integral of ACE as the 

control signal C4-6I1. It has been found that the use of ACE 

as the control signal reduces the frequency and tie-line power 
errors to zero in the steady state, but the transient response is 
not sat i sf actoryC6I] . 

To improve the transient response, linear dcentralised 
Load Frequency Control using Pole-Placement CIST and linear 
optimal control theory C6J has been investigated. The realisation 
of such controllers is difficult, cumbersome and expensive because 
the feedback portion of the above controller is a function of the 
complete state vector of the system. Generally, all the state 
variables are not accessible. Even if state estimation 

techniques are used to estimate the inaccessible state variables, 
the data need to be transfered over long distance. This 



involves additional cost of telemetering. 

The third type of controller called "variable 
s t rue ture— system" ( VSS ) controller has been used by some of the 
investigators C14 & 163. This controller algorithm requires only 
two measurable variables, i.e., frequency deviation and tie-line 
power deviation. This controller is practically as simple as 
that of the conventional controller and can be implemented with 
very little additional cost. However, the VSS controller needs 
the switching strategy which may not be sitTiple for large and 
complex system. To select the best control algorithms the 
requirement rriust be known. This requi rerrients have been discussed 
in the next section. 

l.a REQUI REMEMTS FOR gELECTIWG CONTROL STRATEGY 

The following requirements should be satisfied by the control 
scheme : 

(i) Control loop miust be characterized by a sufficient 
degree of stability. 

(ii) Following a step load change, the frequency error 

should return to zero. This is referred to as isochronous 

control. The magnitude of the transient frequency deviation 
should be minimized, i.e., frequency swing should not exceed 


10.02Hz. 



(iii) The static change in the tie-line power flow following 
a step load change in either area must be zero. 

(iv) The integral of the frequency error should be 
minimized so that accuracy of "synchronous clocks" are 
maintained. The time error should not exceed ± 3 seconds. 

It has been felt that the modelling of each component 
evolved in the control system should be discussed. The modelling 
presented here in brief and is to fimiliarize the reader with 
basic functions involved in the control system. It is given in 
App. A & B. 

1.3 OBJECTIVE 

The need of a fast, effective and economical controller for 
load frequency control has been established in the foregoing 
discussions. The main objective of this is to investigate the 
effectiveness of P.I.D. controller applied to load-frequency 
control of interconnected power system. Also, to compare and 
establish that the P.I.D. controlled power system yields 
relatively much better response than variable structure 
controlled system and approximately similar to responses of pole 
placement technique. However, the implementation of P.I.D. 
controller is much simpler than the above two techniques. 



±. 4 . ORGANISATION 


To meet the objective of investigating the P.I.D. controller 
performance for Load Frequency Control of power systems The 
thesis has been organized in the following manner: 

Chapber Two covers Mathematical analysis of steady state and 
Dynamic Response. 

Chapbor bhrao covers the effect of individual control elements on 
single area as well as two area reheat thermal plant. 

Chapber four covers the effet of proportional, integral and 
derivative control when applied two at a time. The load 
frequency control response has been obtained for two area reheat 
thermal plant. 

Chapber five covers the load frequency control response of two 
area reheat and nonreheat thermal plant and Hydro-thermal plant 
with P.I.D. controller. The load frequency control response with 
VSS controller has been compared. 

Chapber clx covers conclusions and future scope of this work. 



CHAPTER - 2 


MATHEKATICAL ANALYSIS OF STEADY STATE AND DYNAMIC RESPONSE OF 

SINGLE AND TWO AREA POWER SYSTEM^ 

2.1 INTRODUCTION 

Before proceeding into P.I.D. controller system, it is felt 
that the load frequency control of power system with and without 
controller should be analysed for its steady state and dynairiic 
response behaviour. In this chapter the mathematical analysis of 
load frequency control of single and two area power systeins has 
been studied. 

2. 2 RESP<»fSE CtF LOAD FREQUENCY CONTROL OF AN ISCg^ATED 
C SINGLE ARE A3 POWER SYgTEN 

Now we shall present the analysis of load frequency control 
of a simple power system. The corriplete system to be considered 
for the design of the controller may be obtained using equation 
(A. 9), (A. 16) and (A. 26) of App. A. It is shown in figure <E.1). 
The systemi response has been obtained for uncontrolled and 
controlled cases in the following sections! 

2. 2Ca> UNCONTROLLED 

Steady State Response ! In uncontrolled case, speed changer 
has fiKed setting i.e., AP 

^ f f 


ref 


0 




Gh(s) G\{s) Gp(s) 



Figure (2-1) Block diagram representation of single area 
thermal plant (UI'iCa'JTROLLED) 





For step load change 

M 

Laplace transform of it is Pp(s) = -g 
Now from the block diagram (fig. S.1) we obtain by inspection 

^ AfD - APj^> Gp = Af (2.1) 

In Laplace transform form 

Gp 

Af(s) = - :f (2.2) 

Using the Final value theorem 


® ^P M 

Af = lim CS Af(S)> = 7 — x g (2.3) 

s^o 1 +(^)GpG^G^ 

KpM ^ ^ ^ 

“ 1+Kp/R “ ~ D+1/R ^2 
if /? = D + ^ P.U. MWHz 
M 

then ~ ^ constant error (2.4) 

where ft is called area Frequency. Response Characteristics (AFRO. 
Thus in uncontrolled case we find that the steady state response 


has constant error 
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Dynamic Response 

Finding the dynamic response is quite straightforward. By 
taking inverse laplace transform of equation (E.1), gives an 
expression for Af<t). However, as Gj_j, G^, and GpContain at least 
one time constant each, the denominator will be of third order, 
resulting in unwieldly algebra. 

We can simplify the analysis considerably by making the 

reasonable assumption that the action of the speed governor plus 

the turbine generation is instantaneous compared with the rest 

of the power systein, ^t^^ ^p^ where Tp is generally 20 

sec. and T , T*< 1 sec. Thus we can assume T = T.= 0 and gains 
g t g t 

equal to 1. 


Af(s) ~ 


• P 

1+ ^ \ 

^ 1 1+STp 



< 2 . 5 ) 


Equation (2.5) can also be written as 


Af(S) = 


M 


R+Kp S 


S+ 


R+Kj. 

rtT 


( 2 . 6 ) 


Taking inverse Laplace transform of <2. A) we get 


Af (t) 


RKp 

■^*^D R^P 




-e 


t 


r ^P'f^ 

I RTp 



(2.7; 





Thus the error 


persists in uncontrolled case 



A simulated response of uncontrolled single area nonreheat is 
shown in fig. (2.2). The system parameter values are also given. 

2.Z(b) CONTROLLED CASE ; Before analyzing the controlled case of 
a single area power system, we shall define "Control Area". 

CONTROL AREA! The power p>ools in which all the generators are 
assumed to be tightly coupled, i.e. they swing in "unison" with 
change in load or due to speed changer settings. Such an area, 
where all the generators are running coherently, is termed as 
"control area". 

INTEGRAL CONTROL ! By using the control strategy shown in fig. 

we can control the intolerable dynamic frequency changes 
with changes in load and also the synchronous clocks but not 
without error during transient period. We have added to the 
uncontrolled system in fig. 112.13 an integral controller which 
actuates the speed changer by the real power commond signal AP^ 

APc = - K^jAfdt (2.8) 



integral ,APd(s) 

controller 


-Ki 


Kg 




Kp 

S 


H5Tg 




1+STp 


1 

R 


Figure (2-3) Block diagram representation of 
single area (CONTROLLED) power 
plant 




It 


The negative polarity must be chosen so as to cause a positive 
frequency error to give rise to a negative, or "decrease" 
command. 

Here K and K^. are such that K A 1 
g t g t - 

In central load frequency control of a given area, the signal fed 
into the integrator is refered to as area control error (ACE), 

i . e . 

ACE = Af (S.9) 


Taking Laplace transform of equation (2.9). We get 

K. 

AP^(S) = ^ Af(S) 


and for step input load: 


'APjp(S) 


M 

S 


Thus from the block diagram, of fig. (2.3) we get 


K, 


Af(S) ==- 


. K . K 

(l4.STp) + (^ + ^) X 


M 

^ S 



X 


( 2 . - 10 ) 


RK„S( -1+ST )(1+ST. ) 

_ P g t_ 

S( -l+ST. ) ( 1+ST. ) ( 1+STp)R+Kp(RK .+S) 
y ^ r r 1 


M 

S 


using Final value theorem, we readily obtain from the above 
ation the static frequency drop as: 


^ ^ steady 


= limCSAfCS)^ 
S-fO 


0 i.e. no error. 


uiated frequency response of single area non-reheat thermal 
nt is given in fig. <2.4). Integral gain value has been taken 

.;vi C6 :l 


= 0.01 

T = .08 sec 

•3 

= 0.38 sec. 
Tp = 20 sec. 


K.=0.6 

1 

Kp= 120 (Generator Gain) 

K K =1 
g t 

R=2.4 


LIMITATIONS OF INTEGRAL CONTROL 

II has been found C63 that the use of ACE as the control 
nal reduces the frequency and tie-line power errors to 2 ero in 
steady state, but the transient response is not satisfactory. 
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It can be seen in fig. (2.4) that overshoot is more than 0.02Hz 
and setlling time is also more compared to VSS controller and 
proposed controller applied to thesame system, which will be 
illustrated in chapter 6. 

2.4 RESPONSE OF LOAD FRE(5UENCY CONTROL OF TWO AREA INTERCONNECTED 

POWER SYSTEM. 

In this section, analysis of load frequency control of a two 
area power system is presented. The complete system to be 
considered for the design of the controller may be obtained using 
equation (A. 9), (A. 16), (A. 26) of appendix (A) and equation 
<B.8), of appendix (B). 

The system response has been obtained for uncontrolled and 
controlled cases in the following sections. Figure (2.5) 
represents the block diagram of two area uncontrolled 
inter-connected power system. 


UNCONTROLLED CASE : Steady State Response: 

We consider first the uncontrolled case with AP^^= AP^g=0. 
Suppose that the load in each area is suddenly increased by 

Due to the incremental loads, we 
shall have frequency drops in the steady state and these drops 
will be equal to 


incremental steps AP^^ & AP^g. 


Af1 * . = Af2 . . = 

stat stat stat 





Figure (2-5) Bloc 
inte 


?r 



gram representation of two 
inected power plants (Ur>!CCSNiTROLLED) 
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and we will have incremental tie-line power in steady-state, i.e. 


also 


12stat* 

Gistat R. stat 

1 


( 2 . 11 ) 


G2stat 


Rg ^^stat 


( 2 . 12 > 


d(Af1) 


Refering equation (A. 7) for area 1 and setting 


-= 0, we 


write for static conditions 


^*^G1stat ^^D1 " ^1 ^^Istat^ ^^12stat 


(2.13) 


Similarly, for area 2 we can write 


AP__ . AP__ = B_Af_ . ^ + AP_^ . . 

' G2stat D2 2 2stat 21stat 


or AP__ . .- AP__ = B„Af_ . +- AP ^ ^ 

G2stat D2 2 2stat 12stat 


(2. 14) 


since 4P^g s=-APg^ 
or AP^g = a^gAPg^ 


where a. 


Making use of Eqs. (2.10), (2.11) and (2.12) we get the following 


^^Pd2'*~^12^^D1^ 

Z 1 
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AP 


lEstat 


*®2* ^g’ 


Let . B,-. J or (»2 = Bg e 1 

1 H 


Thus we can write 


Af 


stat 




H 


AP 


(/^^APDg-Z^^^P^g) 


Igstat ^S'^^-Ig^ 


P.U.M.W. 


(g. 15) 


(g. 16) 


In case we have B^= Bg= B 



Then (g.15) becomes Af 


stat 


AP +AP 
^^Dg D1 

2fi 


H. 


and (g.16) becomes AP 


^'^Dg ^*^01 


P.U. MW 


Igstat g 

Suppose a step load change occurs in area 1 only, then we get 


Af 


stat 


^'=D1 „ 

Tjr “z 





AP 


AP 


D1 


12stat 


P.U. MW 


Thus we conclude from above that for a load change 
a two area system, the steady state error reduced 
tie-line flow by 50X. 


E . 5 Dynamic Response 
Uncontrol led Case 

Let us consider the transient period. For 
simplicity we shall assume the two areas to 
Further, we shall be neglecting the time constants 
and turbines as these are negligible compared 
constants of generators. 

or T ^ 

P 1 1 1 g 1 

^PE^''" ^tE ^gE 


in area 1, in 
to SOX and 


the sake of 
be identical, 
of governors 
to the time 


For the uncontrolled case AP^^= AP^ 2 = 0, 

following equation by refering to fig. <E.5) 


we 


can write the 


Af (S) 

Afg(S) 


K 


P1 

1+ST 


K 


PE 

1+STP 


Af (S) 

C— i + APp^(S) + 

Af (S) 

C ^ + APjjgCS) - AP^2(S)D 


(E.17) 

(E.18) 



X/ 


also we know CAf ^ (S)-Afg(S) □ 


After substitution we get 


^ ; iT 4»l"Kp 

S +S(-y + ^ )+ L 


2nT K, 


Since ~ 4 HB/P.U.MW Tr, = sec. 

Pi P Bf° 


we gets 


AP^2<S> = 


AP^_(S)-APD,(S) 
D2 1 


o2^/Bf° ^ f° ^ 2nT°f° 
® ■" “H 


,< 0^.0 
nf T 

H 


We see that the denominator is of second order and has the 

„2 „ „ 2 

S + 2oiS + <i>n 


f 1 

where « = ^ 

2 2nT f 
~ H 

2 

If a and both are +ve then system is stable damped 
The roots of the characteristice equation are 


form 



xo 


O 4. / 2 , 2, 

1 ,E ' n 

Now we have three conditions! 

(1) If ot = the system is critically dampled. 

(£) If a > 6)^ the system becomes overdamiped. 

(3) If a < , underdamped. 

A typical ex.a/ri 4 )le „pf two area power systeiTi with paramieter given 
as shown in fig. (2.6) 

CONTROLLED CASE i 

Here also we will use integral control in each area. 

Here in this case there are two control signal variables, in each 
area, namely frequency deviation and tie-line power deviation, 
thus we redefine area control error ACE for the two-areas as 
f o 1 lows ! 


ACE^ = 

^2 


(2.18) 

ACEg= 

*^21 

" V ^2 

(2.19) 


The block diagram representation of two area inter connectedd 
power system is shown in fig. (2.7). 

The speed changer of governor will get the following command 


signals : 























K. . I (AP,_+B.Af . )dt 
i1J 1E 1 1 


< 2 . 20 ) 


AP 


C1 


AP 


C2 



<APg^+B2Afg)dt 


( 2 . 21 ) 


2.6 STATIC RESPONSE 

The choosen strategy will eliminate the steady state 
frequency and tie-line deviations for the following reasons. 
Following a step load change in either area, a new static 
equilibrium, if such an equilibrium exists, can be achieved only 
after the speed changer commands have reached constant values. 
But this evidently requires that both integrands in equation 
(2.20) & (2.21) be zero. 


AP 

AP 


12,0 

21,0 


+ B^Af 
1 

BgAf 


o 

o 


0 

0 


But as we know ^^21 

Therefore Af = AP_, = 0. 

o 21,0 

We conclude that for the above equations to be satisfied, both 
the steady state frequency and tie-line deviations should vanish. 
When this aim is achieved the system frequency equals the desired 
value and the interchange schedule is met which is the guiding 
principle in pool operation. Each area in steady state absorbs 



its own load. Simulated result of previous example with integral 
control is shown in fig. (2.8). 


2.7 CONCLUSIONS 


In this chapter we have tried to give a balance 
for the analysis of load-frequency control of the 
using the methods of modern control theory. 

Under normal operation the controllers keep 
operating around a pre-selected "normal" state 
excursions. The dynamic models of all controllers 


d presentation 
power system 

the generator 
with minimal 
that have been 


discussed have therefore been linear. 



CHAPTER - 3 


STUDY OF PRCtf*ORTI(»tAL INTEGRAL (CONVENTIOMAL) AND DERIVATIVE 
CONTROL ON POWER SYSTEMS, AS USED INDIVIDUALLY 

3.1. INTRODUCTKW 

In this chapter, the load frequency control response of 
single and double area thermal plant will be studied with 
controller used as proportional or integral or derivative 
controller alone. Simulated results have been shown for study- 

3.2. PRC«>C«tTI<»IAL CONTROLLER 

This controller has been evaluated for the specific example 
of load frequency control problem. It has been already 
established that only proportional controller when used alone as 
controller, is not good. It gives steady state error if value 
of gain is increased, the steady state error reduces but 
transient oscillations increase. In fig. I"3.1a3 the frequency 
response of single area reheat thermal plant is shown for Kp= 0.5 
when controller is in forward path. 

The value of system pararrieters have been taken from C14II. 
For ready reference it is given in table (3.1) 





Table 3.1 


T 

= 0.08 

Kd 

= IEOH /P.U.MW 

9 


P 

z 


= 0.3 sec 

K 

r 

= 0.5 


= 10 sec 

R 

= E.4H /P.U.MW 
z 


= 20 sec 


= 0.01 


3. 3. CC»IVEWTI <»IALC I NTEGRALD CC^TTROLLER 

This is the only controller which brings the steady state 
error to zero. It integrates the error and gives control action 
to speed changer mechanism in governor system so as to bring 
output error to zero. This controller have been in use in Load 
frequency control of power plants for quite long time 01,3,4,53. 
The simulated results of load frequency control of single reheat 
thermal plant with integral controller is shown in figure 
<3. 1-b) . 

3. 4 . DERIVATIVE CONTROLLED 

This controller alone is not suitable for control action. 
As it can be seen in fig. (3.1c) it gives constant error. This 
controller action can be explained as follow: 

When the rate of change of error is higher, the controller 
is more active and as the rate of change goes lower, the 








controller action reduces. If error remains constant (i.e. no 
change in error) then the controller output will be zero. 
Therefore, the Load Frequency control developes a constant error 
response . 

3.5 EFFECT ON TWO AREA I NTERC(»INECTED POWER SYSTEM 

3.5.1 Proportional ConlrollT i The effect of proportional 
controller for different values of Kp = 2.5, 1.5 & 0.5 on two 
area reheat thermal plant is studied. The frequency response. 


when controller is 

in forward path. 

is shown in 

fig. (3.2). 

The 

frequency response 

with proportional 

cont rol ler 

in feed back 

path 

is shown in fig. 

(3.5c). The 

value of 

parameters 

of 


interconnected thermal plant is given in table (3.2) . 

Table (3.2) 


T 

T 

T 

R 

B 


P1 ^P2 

Kp^=Kp2=120Hz/P.U.HW 

r^=Tr2=10sec 

K .=K „=0.5 
r 1 r2 I 

t^=Ttg=0.3sec 


^ =Rg =2.4H2/P.U.MW 


. =B„ =0.545 

1 2 

AP^^=0.01 














24 


3. S. a I ntegr al Conlr ol i Integral control in forward path, with 
gain value 0.3, 0.67and 1.0 have been used to study the load 
frequency control response of two area-interconnected power 
plant. The responses have been shown in fig. (3.3). The 
response with = 0.67 in feed back path is shown in fig. 
(3.5b). It can be seen that integral control in feed back path 
is not good, it gives steady state error. A better result is 
obtained with controller in forward path with gain value 0.67, 
which has been also found by Fosa & Elgand C63 . 

3.5.3 Derlvat-lve Conbrolleri t This controller alone in forward 
path or feed back path does not give good response. Simulated 
responses for gain value = 0.05, 0.5 and 1.0 in forward path 
has been obtained. The frequency deviation response is shown in 
figure (3.4). The same with Controller in feed back path with 
gain value 0.5 is shown in fig. {3.5a). 

3.6 C<»ICLUSIONS 

We can conclude from the foregoing results that only 
conventional controller can be used for load frequency control if 
we are required to use any single controller. In fact only this 
controller have been used in the past. 



CHAPTER - 4. 


STUDY OF PROPORTIONAL^ INTEGRAL AND DERIVATIVE CONTROL 
W HEN ANY TWO CC^fTROLLERES ARE USED SIMULTANEOUSLY 

4.1 INTROEHfCTION 

In this chapter, we shall study the effect of two controllers 
used in load frequency control of single as well as double area 
power plants. Two controllers together, normaly proportional 
plus derivative, proportional plus integral and integral plus 
derivative are used in forward path as well as feed back path to 
arrive at appropriate control configuration which yield's best 
response . . 

4.2 PR<W*ORTIOWAL PLUS MIRIVATIVE CONTROLLER 

With this controller used in single area, it is observed 
that it gives steady state error in frequency deviation, refer 
fig. 114.11]. It also gives steady state error in load frequency 
control of two-area interconnected power system. Fig. C4.2all 
shows the frequency deviation when both controllers are in 
forward path. Fig. C4.2b3 shows the response when proportional 
is in forward path and derivative in feedback path. Fig. C4.2cT 
shows the response, when derivative is in forward path and 
proportional controller in feed back path. None of these 
responses are acceptable since they result in steady state error 





Freq.dev.(in hz) 








0.0 



Rg,(4*^ )LF,C. r«sp 
plant with combinatii 
fi9(<3) wth kp.ki in f 
kpin forward ond ki 
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4,. 3 PROPORTIONAL PLUS INTEGRAL CONTROLLER 

Only some controller arrangement gives good response, i.e., 
no steady state error. When proportional and integral both are 
in forward path with Kp = 0.5 and K. = 0.67, the response maximum 
overshoot is 0.024 and sattling time is approximately. 13 sec. 
(fig. 4.3a>. In fig. C4.3bD, when proportional is in forward 
path and integral in feed back path, the maxmimum overshoot is 
O.OES but settling time is approximately 11 sec. When 
proportional is moved to feedback path and integral to forward 
path, both overshoot and sattling time reduces, ( fig. (4.4a)). 
In fig. (4.4b) it can be seen that with Kp in forward path and 
(1+A/S) in both the feedback path, the response has 0.024 
overshoot and settling time 12 sec. When we use integral 
controller in forward path and (1+AS) in both feedback path, the 
response overshoot reduces to 0.012 and transients die out but 
the settling time is approx. 15 sec. as evident from fig. (4.5). 

Therefore, it can be established that I-P arrangement, i.e. 
integral in forward path and proportional in feedback path 
results in the best response. 

Load frequency responses of single area with different 
values of proportional gain and integral gain are shown in fig. 
(4.6) to fig. (4.9). 



Time(sec) 

Ffg.('^-6)_->LrC response of single area thermo! 
plant with integeral+proportional controller 
Ki=0.7 kp=0.3 



0.00 10.00 20.00 

Time(5ec) 

Ffg.(4.7)-+Lrc response of single area thermal 
plant with integerol+proportlonol controller 
Kl=0.3 kp=0.3 





Time(sec) 

Fig.(4‘8)“»LFC response of single area thermal 
plant with inteqeral+proportional controller 
Ki=0.3 kD=1.8 


Fiq.(4.9)-»LFC response of single area thermol 
plarit with Integeral+proportional controller 
Ki=0.5 kp=1.0 
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4.-4 INTEGRAL PLUS DERIVATIVE CONTROLLER 

Since the integral controller brings back the steady state 
error to zero and derivative controller becomes more effective 
during high rate of change of error, this combination of 
controller also gives satisfactory response of load frquency 
control of single as well as two-interconnected power systems. 
Figures (4.10) shows the frequency response of two area reheat 
thermal plant for different configuration of controller. Figure 
(4.11) to fig. (4.14) gives the frequency responses of single 
area reheat thermal plant. 

4.5 C<»ICLUSI(WS 

From the foregoing results, we can conclude that integral 
controller presence in any combination is a must to bring steady 
state error to zero. However the 1-P configuration, i.e., 
integral in forward path and proportional in feedback path yields 


the best response. 



CHAPTER - 5 


P.X.I>._.£PJROPORTIONAL PLUS INTEGRAL PLUS DERIVATIVE > CQiJTPf^ T EP 

APPLIED TO POWE R SYSTTEMS 


5.1 IHTRQDUCTIC»I 

Investigations on load frequency control of power system had 
been in existance since long C6, 14-163. Various methods have 
been suggested, like, optimal cont rol le r , variable-structure 
controller and pole-placement techniques. All these methods have 
certain constraints which has been discussed in chapter-1. To 
overcome these constraints, a P.I.D. controller approach have 
been investigated in this chapter. In 1982 Hiyama C53 has 
suggested P.I.D. controller, but he used it only in forward path. 
In this chapter different P.I.D. controller configurations will 
be studied, which are based on the study and results discussed 
in Appendix 'C ' . 

5. a IMPLEMENT ATI OH 

P.I.D. controller in different configuration applied to 
single area thermal plant, two inter-connected thermal plant and 
hydro-thermal plant is studied. A simple block diagram 
representation of a particular controller configurations has been 
shown in the next section. The value of parameters of power 
systems is given in table (5.1). All other configurations have 
been shown and discussed in appendix *D'. For simplicity a single 







and 


area subsystems, i.e.. Governor (Gg(S)), Turbine (G^(S)) 
Generator <Gp(S)) have been shown as single block (G(S)). One 
important feature of all configurations is that the overall 
characteristic polynomial remains of the same order. 

5.3 EXAMPLE 

One typical example of P.I.D. controller configuration is 
shown in fig. (5.1). In the secondary feed back path, i.e., 
feed back path* a second differentiation and a constant gain 
areused. In forward path only integrator is positioned. With 
this arrangement, the simulated result of the load frequency 
control responses of two area thermal plant are shown in fig. 
(5.2) and fig. (5.3). The response of hydro-thermal plant is 
shown in fig. (5.4). Block daigram representation of a 

typical Hydro-Thermal power systems is given in fig. (5.5). 


TABLE >5.1 


Thermal plant Hydro plant 


Each area rating = 2000 MW 

Base power = 2000 MVA 

R =* 2.4 Hz/P. U. MW 

B = 0.425 P.U. MW/Hz 

T = 0.08 Sec 
9 


Each are rating = 2000 MW 
Base power = 2000 MVA 
R = 2.4 Hz/P.U. MW 
T^ 2 = 0.545 P.U. MW/Hz 
E = 0.425 P.U. MW/Hz 



Freq,dev,(in hz) Freq.dev,(in hz) 



Timersec) 
ise of tv/o 


Fig,(5.2.)L^C response of tv/o area thermal 
plant wth P.I.D controller (ki/s->in 
forward path)ancl{(H-o1st<s) m B1 F.B.Path) 
ki=1.5->kp=2.5,a1=1.5 (kp in'RTF.B.Path) 



00 10.00 

^ Timefsec) 

Rg.(5.3 )LFC response of tv.-6 area thermal 
plant wiih PJ.D controller (kl/s-^in 
forv/ard path)and((1 +a1s*s^ in B1 F.B.Path) 
ki=.67-*kp=2.5,a1 = 1.5 {kp in‘R1 'F.B.Path) 


Freq,dev,(in hz) 



GOVERNOR REHEATER TURBINE I GENERATOR 



thermal power system; 
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T^= 0.3 sec 

T ■= 10 sec 
r 

K = 0.5 
r 

0.545 P.U. MW/Hz 
Tp= 20 sec 
Kp= 120 Hz/P. U. MW 


T = 1 sec., Tr, = 10 sec 
w R 

= 0.5, = 4S-7^.fec 

Tg = 0.513 sec " 

Kp = 120 Hz/P.U. MW 
Tp = 20 sec 


6.4 OBSEStyATlQHS ON RESULTS 

On a close observation of various responses it can be easily 
seen tha^l maximum undershoot or overshoot is less than 0.02 Hz 
and settling time is quiet less. This is as per the 
'requirements of controllor^ discussed in chapter 1. It can be 
clearly observed that the responses with proposed controller is 
better in all respect. All other results of various 
configuration are shown in appendix ’D* . * 

To have a brief observation and study of results, the 
response characteristics have been tabulated and given in table 
(5.2) to (5.4). 

The tie line power deviation, area control error (ACE) and 
frequency deviation of area *2' have also been shown in appendix 
’D’ for Reheat\Nonreheat thermal and Hydro plants. It has been 
compared with Variable-Structure System (VSS) responses. The 
consolidated results are given in table (5. 5). It is also to be 
noted that controller's gain values have been achieved by trail 


and error method. 
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TABLE - S.2 

RESP(»ISE CHARACTERISTICS OF THERMAL SYSTEM 


Conf i 
gura- 

Fig 
No . 

. Maximum 
Undershoot 


Maximum 

Overshoot 

Setl 1 ing 
time in 

t ion 
No. 


Amplitude 

Duration 

Ampl i tude 

Duration 

sec 



in Hz 

in sec. 

in Hz 

in sec 


1 

D.2 

D.3 

0.017 

0.0175 

4.0 

7.0 

0.0025 

0.0 

8.0 

0.0 

13.0 

8.5 

2 

D.9 

D. 10 

0.017 

0.016 

5.0 

1 .2 

0.0015 

0.0015 

7.0 

0.5 

. 15.0 

8.0 

3 

D. 13 

0.0125 

1.5 

0.003 

1.0 

10.0 

4 

D. 16 

0.018 

1.5 

0.0028 

0.8 

11.0 

5 

D. 19 

0.016 

6.0 

0.004 

5.0 

11.0 

6 

D.22 

0.0035 

1.0 

0.0 

0.0 

1.8 

7 

D.25 

0.008 

7.0 

0.002 

9.0 

18.0 

a 

D.28 

0.0045 

5.0 

0.005 

6.0 

11.0 

9 

5.2 

5.3 

0.001 

0.0012 

3.0 

7.0 

0.0 

0.0 

0.0 

0.0 

3.0 

7.0 




TABLE 

-5.3 



RESPONSE CHARACTERISTICS OF NON-REHEAT THERMAL POWER 

SYSTEM 

Conf i 
gura- 

- Fig 
■ No. 

. Maximum 
Undershoot 


Maximum 

Overshoot 

Setl 1 ing 
time in 

tion 
No . 

Amplitude 

Duration 

Ampl itude 

Duration 

sec 



in Hz 

in sec. 

in Hz 

in sec 


6 

D.34b 

0.0012 

5.0 

0.0 

0.0 

5.0 

9 

D .35a 

0.001 

3.0 

0.0 

0.0 

11.5 
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TABLE -5.4, 

RESPCN4SE CHARACTERISTICS OF HYDRO-THERMAL POWER SYSTEM 


Confi- Fig. Maximum 
gura- No. Undershoot 


Maximum 

Overshoot 

Set 1 1 ing 
time in 

lion 





sec 

No. 


Amplitude 

Duration 

Ampl itude 

Duration 



in Hz 

in sec. 

in Hz 

in sec 


1 

D.7 

0.00188 

3.5 

0 . 0035 

5.5 

9.0 

2 

D. 11 

0.00162 

3.2 

0.0035 

7.5 

11.0 

3 

D14 

0.0016 

3.5 

0 . 0055 

10.0 

13.5 

4 

D. 17 

0.0088 

3.5 

0.005 

4.0 

7.5 

5 

D.20 

0.009 

3.8 

0 . 0065 

4.0 

7.8 

6 

D.23 

0.008 

2.0 

0.007 

2.0 

25.0 

7* 

D.26 

0.0013 

3.3 

0 . 0095 

1.8 

13.0 

8 

D.29 

0.0015 

3.3 

0.008 

2.0 

11.5 

9 

5.4 

0.0015 

3.5 

0.008 

2.0 

13.0 

This configuration is 

like case 14 of table 

(C-1>. 





TABLE 5.5 




System Type of system 

responses 

Controller Fig. 
employed No. 

Maximum 
undershoot 
in Hz 

Setlling 
time 
in sec. 

Frequency 

deviation 

Reheat thermal P.I.D. 

V.S.S. 

D.31C 

D.30b 

0.0022 

0.0066 

6.0 

15.0 

area 

2 

Non-reheat 

thermal 

P.I.D. 

V.S.S. 

D.35b 

D.37b 

0.007 

0.006 

8.0 

12.0 



Hydro-thermal P.I.D* 

D31a 

0.014 

8.0 
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Tie line 
power 

Reheat thermal 

P.I.D. 

V.S.S. 

D.31a 

D.30c 

0.004 

0.006 

6.0 

15.0 

deviation 

Non-reheat 

P.I.D. 

D.35c 

0.003 

7.0 


thermal 

V.S.S. 

D.33c 

0.005 

12.0 


Hydro-thermal 

“P'.l.D. 

D.37a 


576 

ACE 

Reheat thermal 

P.I.D. 

V.S.S. 

D.31b 

D.30d 

0.001 

0.018 

4.0 

15.0 


iMon-reheat 

>.I.D. 

D.35d 

0.062 

8.0 


thermal 

V.S.S. 

D.33d 

0.012 

10.0 


5.5 C<»ICLUSION 

The effectiveness of P.I.D. controller for load frequency 
control of interconnected power systems has been studied in this 
chapter. The simulation results (shown in App. 'D* > and tables 
demonstrate the efficiency of the proposed design. It is also 
shown that the proposed controller satisfies all requirements 
stated in chapter 1 . 



CHAPTER 


6 


CONCLUSIONS 


6. 1 GENERAL 

This report has been devoted to the Load Frequency Control 
of single i^rea thermal and hydro plant as well as two area thermal 
and hy d r 0 “ t he r ma 1 power systems, F-or this purp<ose p'roportional , 
conventional (Integral) and derivative control configurations 
in various combinations are used. The systems have been 
simulated on computer and analysed 'EMPIRICALLY*. The following 
conclusion may be drawn: 

(i) The proportional or derivative controller alone does not 
give good response, i.e., gives steady state error in single 
area as well as two area power systems. 

(ii) Only conventional controller gives steady state zero error. 

(iii) When any two controllers are used at a time only integral 
plus proportional and integral plus derivative controllers yield 
zero error steady state response in single area. When used in 
two area power system, these combinations yield good result only 
in some configurations. 

(iv) P.I.D. controller used in single area thermal or 
hydro plant, shows very satisfactory responses. 

(v) P.I.D. controller used in load frequency control of two 
area interconnected power system shows that only certain control 
configuration yields good results. 

(vi) Load frequency control of interconnected power system when 
P.I.D. controller is implemented, yields 1 ow ove r shoo t ( less than 
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0.02H^), reduced transient oscillations and fast 
settling time. This is within requirement of standard 
controller . 

(vii) The proposed P.l.D. configuration can be usedd in 
multi— area (i.e. more than two) interconnected power 
system C6It. 

6.2 FUTURE RESEARCH 

<a> In this thisis results have been obtained 
empirically. The mathematical analysis can be carried 
out on these systems and a formula or theory can be 
estabal ished. This formula might be used for the 
control of any complex plant. 

(b) The P.l.D. controller implementation can be 
investigated on other area of power generation, 
transmission and distribution system. 

(c) Further investigation can be carried out with 
consideration of Generation-Rate-Constraints (GRC) in 
load frequency control of interconnected thermal power 
plants . 

(d) Load frequency control of interconnected 
hydro plants can be investigated with P.l.D. 


control ler . 
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APPEMDIX A 


MCM:>ELLING FC«t load-pkechjency control study 

The load frequency control system consists of various 
subsystems mentioned below. 

<i) Governor system 

(ii) Turbine system 

(iii> Generator and Load system 

A- 1 Governor System 

The governor is the primary tool for the load frequency 
control. Governing system is basically of two types: 

(i) Mechanical-Hydraulic governor 

(ii) Electro-Hydraulic governor/P . I .D . governor. 

A-2 Mechani cal-Hvdr a ul i c Governor 

A schematic arrangement of the main features of a speed 
governing system of the kind used on steam turbine is shown in 
fig. <A.1). 

The main components of governor system are listed below 

(a) Centrifugal Governor : This is purely mechanical speed 

sensitive device coupled directly to the prime— mover. As the 

speed increases the fly balls moves upwards and the point 2 on 
linkage mechanism moves upwords. 
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- *- ^niEJLLfiiJli. it consists of pilot valve (v) and oil 
servomotor (0). A low power pilot valve movement is converted 
into high power level piston valve ( servomotor ) movement . This 
is necessary for closing and opening the steam valve against high 
pressure steam. 

(c) Linkage Mechanism: 1-2-3 is a rigid link pivoted at point 2. 
For point 3—4—5, another rigid link pivot is provided at point 4. 
The function of linkage mechanism is to control the steam valve 
'V . It also provides a feedback from the steam valve movement. 

(d) S peed Changer : It provides a steady state power output 
setting for the turbine. Its downward movement opens the upper 
pilot valve so that more steam can be admitted to turbine under 
steady condition. 

To develop the mathematical representation of the system, let 
us assume that the system is operating under steady-state 
condition i.e. the linkage mechanism is stationary, pilot valve 
is closed, the steam valve is opened by a definite magnitude, 
turbine is running at constant speed and the turbine output is 
balancing generator load. 

Thus the nominal conditions are; 

Power delivered = 

^ j 0 

Speed * w 

o 

frequency “ f 

Primemover valve position = 



Let us change the speed changer to cornmand a power increase 
) . The speed changer iriovement gives rise to the linkage 
point '1’ moves upwards a small distance (AX^). We may 
establish! 

AX. = KAP 
1 c 

The link point "3" will move downward because of Linkage (1-2-3) 

* 

action. Lot it be (AX^)- Further, the link point **4" moves 

« 

downward by an amount (AX^) which moves the piston in pilot 
servo (v), resulting in high pressure oil flow in the bottom of 

piston of servomotor 'O'. The piston moves upwards by an amount 

(AX^) and steam valve opening increases. Thereby increasing the 
torque developed by turbine. This increased torque increases 
the speed of generators i.e. frequency (Af>. 

This change of speed results in the outward movement of flyball 
of the governor. Thus the link point '2' moves slightly upward a 
small distance (AXg). 

Due to the movement of link point 2, the link point 3 also 

moves upward by an amount AX^ which is also proportional to Af . 

Thus the net movement of link point 3 is 

AX„ = AXA + AX" (A.1) 

3 3 3 

where 

^^3 = - 

= - K. APc 
1 


o r 
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and 


AX^ = Kg Af 


Thus we can write: AX_ = -K.APc + K„Af 


Again 


(A. 2) 


AX. = AX! + AXV 
4 4 4 


<A.3) 


where AX^Clg^+l^^) = AX 2 ( 1 ^^) 


-^^4 ^^ 34 ‘^^ 45 ^ ^^ 5 ^ ^ 34 ^ 


Thus we can write (A. 3) as 


AX, = K3AX3 ^ K,AX5 


(A. 4) 


Now if we make our assumption that the flow of oil into the 
servo-motor is proportional to position AX^ of the pilot valve 
V, then the movement AX^ of the piston (or valve opening) can be 
expressed as 


AX,, = AX 
5 V 



(-AX^)dt 


(A. 5) 


Taking Laplace transform of eqn. (A.E), (A. 4) and (A. 5) we 


get 



Gg (S) 



Figure f/\-2) Block diagram representation of 
' governor 


I 

AX 5 (S) 


G^ (5) 



Kt 

^ 

— — ^ 

1 + STt 



APt (S) = APq{S) 


Figure (A-7) 


.ock diagram representation of 
linearised model of nonhxeheat 
irbine 





AX„(6) = -K APc(s) + K^AF(s) 


(A. 6) 


AX^(s> = K^AXgCs) + K^AX^(s) 


(A. 7) 


Xj^(s) = 


-K- - X. (s) 
5 s 4 


(A.S) 


eliminating AX 2 <s) and AX^(s), we obtain the following equation 


AX^(s) 

3 


K^KjAP^CsI-KgK^AFCs) 

<'<4 " > 

s 


or 


. K 

AX„<s) = CAP (s)-'^ AF(s)3 x 

5 c R 1+sT 


(A. 9) 


V 

where R = jr- = speed regulation of the governor 
2 


K K 

A ^1*^3 

K = —n — = gain of speed governor 
9 n . 


T = = time constant of speed governor, 

9 K4K3 


In standard block diagsram symbol we can represent equation (9) 
as in fig. (A. 2) . 

A-3 Modified Models of Mechanical Governors 

Many authors have considered various other parametersof 
generating station to consider different requirements. For 


example Chowdhary C103, Hovey C123 have modelled for stability 
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(A. 10) 


(A. 11 ) 


where 6 is temporary speed droop 
O' is permanent speed droop 
z is incremental gate position 
is dashpot reset time 
w is incremental machine speed. 

H.E. Coles C113 modelled steam terboalternator by linearising the 
nonlinear differential equations. The transfer function relating 
throttle-valve opening to changes in set speed has been obtained 
by linearising these equations giving 





(1-H- P) ( 1+T .P) ( 1+T_P) 

y 1 2 


where = governor relay time constant 
= primary relay time constant 
Tg s= secondary valve time constant 


(A. 12) 




Figure (A#3) Mechanical-hydro governor for 
steam turbine 



Transient d 
compensation 


Figure (A-4) Mechanical-hydro governor for 
steam turbine 
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constants which relates the steam valve lift to the 
change in speed, at any perticular operating load. 
z = throttle value displacement 
= machine speed rad/sec. 



In 1973 IEEE committee report was published on dynamic models for 
steam and hydro turbines. As per this report steam turbine 
mechani cal-hydroaul i c governor was represented as in fig. (A. 3) 
and ( A . 4 ) . 

A. 4 ELECTRO-HYDRAULIC SPEED GOVERNOR 

An e 1 e c t ro— hydraul i c speed control mehanism provides 
flexibility through the use of electronic circuits in place of 
mechanical components in the low power portions. 

A typical Electro-hydrauti c is PID governor represented as 
in fig. (A.5) . 

Hagihera et al CRef.D have modeled the governor with other 
systema with following assumptions! See fig. (A. 6). 

(1) The hydrosystem operates in linear mode with small load 

perturbations . 

(E) The water column is inelastic 
(3) the load is isolated.. 

where = governor permanant speed droop 

S = temporary speed droop 

T = water starting time 
w 




Permanent 

droop 


Figure (a-^)) Electro-hydraulic governor 
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Figure (a- 6) Electro-hydraulic governor 
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~ mechanical starting time 
a = load self regulation factor 

= dashpot time constant of temporary droop governor. 

If we adjust the gains Kp , H. and of the PID governor to the 
values as below: 

^i= 6T“ ^d " ^ 

r 

Murthy & Hariharan C83 has shown the Temporary droop governor 

i s s 

% = S ■ ' 2T • Xd = ^ 

r 

where is derivative feed back in temporary droop governor. 

A. 5 MODELLING OF TURBINE 

STEAM TURBINE : All compound steam turbine s.y£.tcms . u.tilize 
governor-controlled valves at the inlet to the high presure (or 
very high pressure) turbine to control steam flow. The steam 
chest and inlet piping to the first turbine cylinder and reheater 
and cross over piping downstream all introduce delays between 
valve movement and change in steam flow. The principal objective 
in modelling the steam system for stability studies is to account 
for these delays. 

The continuity equation is 

4 ^ = Q. - Q t 

in out 


(A. 13) 
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where W is the weight of steam in volume V(ft^), t is time in 
seconds and Q.^ and are flows assuming the weight flow out 

of the vessel is proportional to pressure in the vessel. 


Q ♦ ”* *5 ^ 

out P o 


out o dP 

”31“ " dTt 

o 

where P is variable vessel pressure and P^ is steady state vessel 
pressure, is steady state weight flow out at P^ 


dW _ dP _ d dP 


(A. 14) 


where V is specific volume (ft /wt) of steam in vessel 


P dQ ^ 

a. - a = _ V (^. (i> 
in out Q <»P V dt 


(A. 15) 


T = £ V (— ) 

Q ^ dp V' 
o 


Then 


Q. - Q . 
in out 


out _ 1 

15. ” 1+ST 


(A. 16) 


This model is without reheater and shown in figure (A. 7) 
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H.E Coles C11D has modelled steam turbine as 


steam pressure AP* G 

^ 

throttle-valve opening Az "" T+t'T 

where is intrained-steam gain constant. 



where = steam pressure at the inlet to the high pressure 
2 

turbine, Ib/in 

Gg is Reheater gain constant 
z = throttle-value displacement, 

T = entrained-steam time-constant 

T 

T = reheater time constant 

r 

T^= primemover torque P.U. 



Both (A, 16) and (A. 17) are same except the parameters of Reheater 
is included in (A. 17). If there is no radiation, i.e.,n=0 then 
(A. 16) becomes equal to (A. 17). 
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Mode 1 1 in q of H ydr o-turbine 

The transient characteristics of hydroturbine are determined 
by the dynamics of water flow in the penstock. The conversion of 
flow and head (h) to power by the turbine involves only 
non— dynamic relationships. The most precise models of water 
pressure and flow in the penstock are those which treat the 
travelling-wave phenomena. But usually it is not necessary tf 
use travelling wave models for stability studies, although they 
are used regularly for detailed plant design studies. 

The approximate linear model desired by Murthy and Hariharan 
CST & IEEE committee report C9I1 , Hoavey C103 and Chowdhary C1SI1 
is given below: 

The turbine flow is given fay 

q^= g + 0.5h (A. 18) 

and torque M = g+1.5h (A. 19) 

where g is gate position and h is head of water. 

For inelastic water column the head and flow are related by the 
impedance function, 

— = T S where T is water starting time 
q^ w w 

by algebraic manipulation we get 

M 1-T S 

M _ w (A. 20) 

g " TTrrT^ys 

^ w 


is the turbine transfer function with long penstock and elastic 
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water column. 

Murty et al C8D have derived modified transfer function 


„ 1-(T /T )tanh(T S) 

M we e 

g * -| + <T /ST ytanh<T 5 ) 

^ we e 


<A.21) 


where pipe line reflection time. 

Again C.K. Sanathanam C71 have considered the detailed 

Hydroturbine as: 

^ 1-1.9865+0.36065^ 

I SZ ..I. 

1+0.99325+0.36065^ 


model of 


<A.22) 


A. 6 GENERATOR LOAD MODELLING 

We shall develop mathematical modelling of an isolated 
generator, which is only supplying local load and is not 
supplying power to another area via a tie-line. 

Suppose there is real load change of AP^. 

Due to the action of turbine controllers, the generator 
increases its output by the amount AP^. The net surplus power 

AP - AP will be absorbed by the system in two ways. 

G D 

1. By increasing the kinetic energy W^.^ in the rotor generator 

at the rate d(W . ) 

C 

dt 

2. By an increased load consumption. All typical loads 
experiences on increase 

B « ygH MW/H^ 



51 


B is called damping coefficient. 

Thus we can express this surplus power 


AP^-AP 

u 





Kin 


+ BAf MW 


<A.E3) 


As the kinetic energy is proportional to the square of the speed, 
we can thus write 


W 


Kin 




M.W.S 


as f = 
frequency 


f o+A f 
after 


where f° nominal 
disturbance. 


frequency. 


and f is new 


W 


Kin 


*^Kin' ,0 


W 


Kin 


C1 + 


2Af 

„o 




^ ( 1+2 
Kin 


‘^^‘^kin^ ^^Kin dCAf) 


dt 


IT 


Substituting (A. 24) in (A. 23) we get 


AP_-AP_= + BAf MW 

G D dt 


(A. 24) 


(A. 25) 


By dividing this equation by the generator rating and by 

introducing the perunit inertia constant 




Figure (A-8) Representation of generator 
load model 




Speed governing APqCS) 



Figure (A-'9) Block diagram representation 
of power system 
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H 




Kin 


■p— 

r 


MWS/Mw 


(A. 25) becomes 


AP _ AP - 2H dAf 

dT" P®*' MW 


Taking the Laplace transform. 


pu 

APq(S)--APjj(S) = ~ SAf(S) + BAf(S) 


or we can write Af (S) = G„(S)CAP^(S)-AP„(S)D 

P G D 


(A. 6) 


A ^P 
= TTST 


— — - ^ 

~ f°B 


sec . 


Tp = ]g H^/P.U. MW 


Thus the block diagram representation will be as in fig. <A.8). 

A. 7 Block Daiqram Representation of single area power system 

We can easily obtain its model by combining the block 
diagrams. Shown in figure (A. 2), (A. 6) and <A.7). By properly 
connecting these block diagram, we get fig. (A. 9) which is a 
single area representat ion . 
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APPENDIX jg 

!l Q ] D.iLXMG . OF INTERCONNECTED POWER POOLS 


B.1 LOAD FREQUENCY CONTROL OF TWO-AREA SYSTEM 

FrofTi ^ prscLic^l poinL oP vi©Wy Ih© prol)l6Fri PrGQUfincy 

control oT inter~connected sreais, or power pools, are niore 
important than those of isolated areas. Practically, all powpcr 

systems today are tied together with neighbouriq areas, and the -rS 

...... 

problem of load-frequency control becomes a joint undertaking. 
Closely associated is the problerri of controling the power flows 
on the inter-ties. 

Since peak demands occurs at various hours of the day in 
various areas, the ratio between peak and average load for a 
large pool is smaller than that of the individual systems. It is 
obvious, therefore, that all pool members can benefit from a 
reduced need of reserve capacity by a scheduled arrangement of 
energy interchange. 



where 6° and <5° 
1 2 


are 


the angles of end voltages and Vg 
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respectively. The order of the subscripts indicates that the 
tie-line power is defined positive in direction 1 to2. 

For small deviation in the angles 6^and 6^ the tie-line power 
changes with the amount 


APig % cos (6^ - 6") (A<5^-AAg) 


(B.2) 


Analogous to the concept of "electric stiffness" of 
synchronous machines we define the "synchronizing coefficient" 


of a line 


o A I 1 N '^2 1 o o, 

T = cos(A®-T5p MW/rad 


(B.3) 


Thus we can write equation (B.2), tie-line power deviation 


AP,„ = T (A<5 -A6^) MW 
12 12 


The frequency deviation Af is related to the reference angle AS 


by the formula 


1 d .rO 

m ' 


1 (AS) 
2n dt 


d(Ai5) 


=* HnrAf 


A45 = 2n I Af dt 
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respectively. The order of the subscripts indicates that the 
tie-line power is defined positive in direction 1 to2. 

For small deviation in the angles 6^and the tie-line power 
changes with the amount 


AP 


12 




y -O -O . 

cos i6 ^ - 6^) 


(AA^-AAg) 


(B.2) 


Analogous to the concept of "electric stiffness" of 
synchronous machines we define the "synchronizing coefficient" 
of a line 


T “ cos(A^-ig) 


MW/rad 


<B.3) 


Thus we can write equation {B.2), tie-line power deviation 


AP^2 "" T°(A<5.^-A6g) MW 

The frequency deviation Af is related to the reference angle AA 

by the formula 

i H o ^ (A<5 ) 

^ ^ (S +AS) « dt~ 


or 


d(A6> 


a EfiAf 


or 


A^ « 2fi 



0 


Af dt 




Figure (B-1) Representation of tie line 
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By expressing tie-line power deviation in terms of Af rather than 
A<5 we thus get 


APi2 


t 

&tT°(J Af^dt 


t 

J Afgdt) MW 


<B.4) 


Laplace transformation of the last formula yields 


A P ^ _ ( S ) 
12 


2nT 


(Af^ (S)-Afg(S) ) 


(B.5) 


Representing this equation in terms of block diagram symbols 
yields the diagram is given in fig. (B.1). 

Similar like equation (5) the incremental tie— line power expected 
from area 2 is given by 


AP„,(S) « 
21 


2nT® 


<Af2(S)-Af^ (S) ) 


(B.6) 


The power balance equation for single area is 


AP - AP 
G D 



d<Af^> 


+ BAf^ 


Now in case of two area we must modify the same equation for area 
1 and area 2 accordingly. 

Thus we must clearly add right side of this equation 

that is. 


AP_ .“AP„ 
G1 D1 


dCAf^) 


"ar 


^ B^Af^^ AP^g 


(B.7) 


Taking Laplace transform of this equation 




B-2) Block diagram representation of two inter-connected 
pov.C'i plant 










where GP^(S) = l^Tp^ ‘ 

Thus the complete Block diagram representation of two area LFC is 
as given in figure (B.2). 
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appendix c 


OF A FIRST ORDER LEAD-LAG 

SYSTEM WITH P. I.D. CPRCW^ORTKWIAL PLUS INTEGRAL PLUS 

DERIVATIVE5CONTROL 


C.l C^IMERAL 

In literature , the performance and design procedures for 
P.I.D. controller have been discussed C3 & 51 . However, various 
methods to implement P.I.D. controllers have not been given much 
attention. In this appendix, an attempt has been made to 
determine the possible methods for implementing the P.I.D. 
controllers, and selecting the feasible methods. The selection 
of feasible methods has been done by determining the steady— state 
error. Normally, it is desired that the steady-state error must 
be zero for given control scheme. To investigate the feasible 
methods for implementing the P.I.D. controllers, the system 
transfer function considered is given belowi 


G(S ) »= 


1-T S 
w 

ITT 'g 

w 


The criterion for choosing the configuration for the PID 
scheme has been based on characteristic equation, i.e. same 
characteristic equation must be obtained for various 
configuration, that is: 

P(S *= 1+C<K„+K .S+K. /S)G(S)3 
p d i 


















*^o 


C.2 SYSTEM REPRESENTATION 

Fig.<C.1) shows the general block diagsram of the system. 
The block G^^<S), G^^(S), 6^3(8) and G^^(S) is the various 

controller implementation positions. In any one configuration, a 
single control block or combination of blocks can be used. 
Further, a single block may contain single control element or 
combination of three controller, proportional, integeral and 
der i vat i ve « 

For illustration, an example is shown in figure (C.E). In 

this configuration 

G^^<S) = Kp+K./S+K^S and 

G „(S) = G „(S) « G „(S) = 0 
c2 c3 c4 

To analyse the steady state error of this system we can proceed 
as follo-ws* 

It can be found that 




From Final value theorem, when R(S) = ^ (step input) 

(1) steady state error >= lin, SE<S) 

S^O 


1 im 


&*0 1+(Kp+Kj/S+KjS)G(S) ® 
0 




= 0 


(C.2) 


( 2 ) 


Steady state error with R(S) ~ 1 (Impulse input) 


* lim SE(S) 

S,0 


, . S 0 

»= 1 1 m • — — — =-<- — = 0 

S-0 1+(Ktj+K,/S+K,S) ^i 

r 1 d 

C.3 BRIEF STUDY 

It can he seen from the foregoing example that likewise 
number of configuration can be studied. To avoid exhaustive 
study, the configuration and results has been presented in 
tabulated form. Refer Table (C.1). Total 28 configuration were 
investigated and 14 were found satisfactory i.e. steady state 
error zero. It is to be remembered that in all systems, the 


characster isti cs polynomial remains same. 
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TABLE C.l 



No. when when 


R(S)=1 R<S)= 1 

s 


1 . 

K„+K . /S+KdS 

P 1 

Open 

1 

open 

0 

0 

E. 

K. /S+KdS 

1 

^P 

1 

open 

0 

0 

3. 

K./S 

X 

Kp+KdS 

1 

open 

0 

0 

4. 

K . /S+K„ 

1 P 

KdS 

1 

open 

0 

0 

5. 

K./S 

1 

open 

1+As+AgS^ 

open 

0 

0 

6. 

Kp 

open 

A ^4*AgS 

KdS 

0 

0 

7. 

K . /S+K 

1 e 

open 

1+AS 

open 

0 

0 

a. 

K./S 

1 

KdS 

1+AS 

open 

0 

0 

9, 

K./S 

1 

Kp 

1+AS^ 

open 

0 

0 

10 . 

K,/S 

i 

open 

1+AS^ 

•<P 

0 

0 

11. 

K./S 

1 

open 

1+AS 

KdS 

0 

0 

IE. 

1 

open 

KdS 

KpfK./S 

0 

0 

13. 

KdS 

0 pen 

A .+A„/S 

1 S 

K./S 

i 

0 

0 

14. 

Kp 

open 

1+AS 

K./S 

1 

0 

0 

15. 

1 

open 

K./S 

1 

Kp+KdS 

0 

1 

16. 

Kp K 

. /S-t-KdS 
1 

1 

open 

0 

1 

17. 

Kp+KdS 

K./S 

1 

1 

open 

0 

1 

18. 

1 

Kp 

K./S 

KdS 

0 

1 

19. 

KdS K 

./S+Kp 

1 

open 

0 

1 

SO. 

KdS 

open 

Ai A3 

open 

0 

1 

SI. 

^P 

open 

open 

0 

1 



0 


1 


A. A, 

^ . /V . ’ 


21. 

Kp 

open 


open 

22. 

KdS-+Kp 

open 

■g— » 

open 

23. 

KdS 

K. /S 
1 


open 

24. 

1 

KdS+K| 

K , 

P F" 

open 

25. 

. 1 

A^S 

open 

26 . 

open 

d P i ^ 

K . /S 

i 

open 

27. 

K .S 
d 

open 

A^+Ag/S 

Kp 

28. 

1 

KdS 

K^/S 

Kp 


C.4 OBSERVATIONS 

From the study of above all configurations, we find that 
only 2fi configurations are possible, keeping the order of 
characteristic polynomial same. Out of these possible 
configurations, 14 are found to be suitable to yield steady state 
error zero with step input. 






APPEWDDIX - ^ 

C<»miOL CONFIGURATIONS 

The various controller configurations are illustrtated 

belowi 

D.1 CONFIGURATION ONE 

The controller is simply positioned in forward path, as 
shown in fig, (D« 1). The frequency response of two area reheat 
thermal plant with P.I.D. controller is shown in fig. (D.E) and 
(D.3) for different value of gains. It can be compared with 
frequency response with conventional controller in fig. (D.4). 
In fig- (D.5) and (D.6), the freuency response of single area 
thermal plant is shown with P.I.D. controller and conventional 
controller respectively. The response of hydro-thermal is shown 
in Fig. <D.7). 

D.a CONFIGURATION TlfO 

For this case the P.I.D. controller configuration for the 
power plant is shown in fig. (D.fi). 

The load frequency control response of the power system in the 
above configuration, is shown in fig. (D.9), (D.10) and (D.11) 
for different value of gains. 

D.3 CONFIGURATION THREE 

In this case the P.I.D. components have been arranged as 
shown in fig. {D,1E>. The load frequency control response of 
this system i« shown in fig- (D.13) and (D.14). 














Freq.dev (area 1) Freq.dev (area 1) 





rre<iMmX 



2 

Wfi'i 


D.02C' 


0 OtK i 


0.020 ■■ 


0,040 


- 0 . 0 > f)0 



I j |-'r-r”f"rT'T“rT ,-,-7-T-rT~rT’i 

10.0 20.0 30.0 

• Tirrir (f.qr) 


of tv^o or«n tli«rn'iDl 
plont vvi^h {'*Ofl7EJif'IONAL controller ki^O.67 



Freq.dev.Cin hz) Freq.devcm i 



0,00 10.00 20.00 


Tinne(sec) 

'y CD -5, ) lFC Respono# of Single Area Thermal Plant 
0*^P1D control b- conventional (int.)control 



Time(sec) 

Fiq.(D-6)LFC response of one area reheat->thermal 
plant with PID controller l<p=1.8.ki=0.5 and kd-0.5 



Freq.tfev,(in hz) 















0.01 



-o.o:$ 


0.00 


r I I , 1 I |-T-r-T-_ 

„ 20,00 30X0 

. , . . Time (sec) 

ig ^(P' ' ^ ^ of' two area roheat-^thermol 

plant y*‘t(' MD controller 'Iki+kd fn forward path) 
kp-b.O ki“0.5^ ond kd=5.04(kp fn feedback poth) 



- 0,03 - ■ 

1 ^ f 7 
0,00 


■TT" I I " I [ ' r -i'-i I I I I I I I I I I r I 

10.00 20,00 30.00 


Time(sec) 

ng.ft ),10 )LFC reeponsfl of two area reheot-+thermal 
plant with PID controller (ki+kd in forward poth) 
kp=1.6 ki=0.67 and kd=0.05(kp in feedback path) 



eq.dev.( 



P'9 fB-H) (-♦LFC response of two area 
plont with P.I.D controller (ki/s+kd 
forward path )and (kp in T.b.path) 
kd^O.B-^Kp-O.S ki*0.57 


THERMAL HYDRO 
s in 



Freq.ciev.(ini hz) Freq.dev (areo 1) 




0.01 


0.00 


0.01 


~Q.Q2 







D.4 CONFIGURATION FOUR 

The «rr«»ng#nrenl of P.l.D. controller is shown in fig. (d. 15 ) 
and the load frequoTicy control response is shown in fig. (D 16) 
and (D.17). 

d.b configuration five 

T)»is is tl.e it*ost compieK of all the systems due to the 
double di f f er ent iat 1 ott in the secondary feedback path. The 
implementat J on i ». sl.owjt in fig. (D. 18). The load frequency control 
response is shown in fig. <D.19) and (D.20). 

D.6 CONFIGURATION STX 

In this tact* both the feedback path and forward path have 
been impl en.ented with components of P.I.D. controller as shown in 
figure (D.?1). Tltr load frequency control response is shown in 
fig. (D.22) and <D.23>. 

D.7 CONFIGURATION SKVEN 

In this controller configuration, a differentiation is 
positioned in secondary feedback path and effective P.I.D. 
controller is i rtipl emented as shown in fig. (D.24) and the load 
frequency ri‘sputt‘*e for different value of integral gain, is shown 
in fig. (D.25> aird 
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D.e CONFIGURATION EIGHT 

This conf Kinration , P.l.D. controlller is shown as in fig. 
(D.27) load frequcMicy control response is shown in figure 

(D.28) and (0.29) for different Ki & Kd. 

D.e CONFIGURATION MINE 

It is shown and discussed as an example in chapter 5. 
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